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Abstract 

It is shown how strictly four-dimensional integration by parts combined 
with differential renormalization and its infrared analogue can be applied 
for calculation of Feynman diagrams. 
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1. Integration by parts (IBP) [[j] within dimensional regularization ^] turned out 
to be one of the most powerful methods of calculation of Feynman diagrams. In spite 
of the fact that dimensional regularization and renormalization || are commonly used 
in practical calculations they happen to be inconvenient in situations where chiral and 
super symmetries are involved. The purpose of this brief note is to present a strictly 
four-dimensional version of this method of calculation. It will be shown how (usual) 
four-dimensional IBP supplied with technique of differential renormalization [4-9] and 
its infrared (IR) analogue (-R-operation) can be used for calculation of diagrams. 

In the next section we shall describe prescriptions of differential renormalization and 
i?-operation for simplest diagrams. Using an example of the master two-loop diagram, 
we shall show in Section 3 how four-dimensional IBP is applied for calculations. In 
conclusion possible applications of the proposed method are applied. 

2. Differential renormalization was first defined in coordinate space [4-8]. It is pos- 
sible to translate its recipes into the momentum space language. Another possibility 
is to use homogeneity properties of Feynman integrals and formulate the corresponding 
recipe just in momentum space — see where this prescription was formulated for loga- 
rithmically divergent diagrams with the simplest topology of subdivergences. (Probably, 
it is more natural to call this renormalization homogeneous rather than differential.) 
In particular, for the simplest one-loop propagator-type Feynman integral^ one has the 
following renormalization prescription J9| : 

Here /i is a massive parameter that takes into account finite arbitrariness of renormaliza- 
tion. The derivative in q is applied before integration so that the corresponding integral 
happens to be convergent. 

To calculate (|j) it is worthwhile to introduce analytic regularization ]lOj : 



RF(q) = / d 4 fc -— -q 



Id [i 



2 A 



dXJ 2 H dq{k 2 ) 1+x {q-k) 2 



(2) 

A=0 



When A ^ 0, we may use the following order: to calculate the integral, differentiate in 
q, differentiate in A and finally put A = 0. When calculating the integral one uses the 
four-dimensional one-loop formula 

d4fc (j- fc)2)A 2 = ^ G ^ X ^ ((? 2 ) A 1 1 +A 2 -2 > ( 3 ) 

where G is the four-dimensional G-function 

rn ^ ,_ r(A 1 + A 2 -2)r(2-A 1 )r(2-A 2 ) 

G(Al ' A2) " r(Ai)r(A 2 ) f(4-a 1 -a 2 ) • (4) 



For simplicity, we consider Feynman integrals in Euclidean space. 
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In particular, 

G(M + A) = ^. (5) 

Finally we have 

R J dik k^ L w =n ^ 1 - lnq2/ ^- (6) 

We shall need as well an IR analogue of differential renormalization. For the case of 
dimensional renormalization, the corresponding IR analogue was developed in This 
operation R removes IR divergences (off the mass shell) in a quite similar way as usual 
dimensional renormalization removes ultraviolet (UV) divergences. When we combine it 
with the i?-operation itself we obtain the so-called i?*-operation Jll| (R* — RR) which is 
successfully applied in renormalization group calculations and in asymptotic expansions 
of Feynman integrals in various limits of momenta and masses (see. e.g., for a 
review) . 

Differential i?-operation can be defined in momentum space in quite the same way as 
differential renormalization in coordinate space.^] For the simplest IR-divergent subgraph 
consisting of two subsequent massless lines that generate a nonintegrable factor 1/k , we 
have the following prescription: 

where p. — is an IR-renormalization parameter. Therefore, for the one-loop IR-divergent 
Feynman integral 

d 4 fc, ' 



'fc 4 (g- fc) 2 

we obtain the following 'i?-normalized' expression: 



R / d 4 fc— — = / d 4 fc -—k — - — . (8) 

J k 4 {q~k) 2 J \2dk k 4 J {q-k) 2 W 

The derivative involved in (]?]) and is understood in the distributional sense (that is 
formally equivalent to the four-dimensional IBP), i.e. 

To calculate (||) it is helpful to introduce analytic regularization and then differentiate 
the integrand at A ^ using 2M fc (fc^) 2 + A = ~^ (>S) 2 +>- • We ^ ave 

k S dik i^w=^ {l+wr " 2) - (10) 



3 To be more precise, renormalization prescriptions in coordinate space formulated in refs. 
reformulated as 'IR-renormalization' prescriptions in momentum space, using UV-IR analogy (as it was 



done in the case of dimensional fi-operation — see jnj [l^ ]) 
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We shall need also similar integrals with additional logarithms which are as well 
calculated with the help of analytic regularization and one- loop integration formula (||). 
One gets 



Ail \n{{q-kf/^) f 4 (I d , Mky^) \ ln( g -fc) 2 / M 2 
' ' d * fc%-fc) 2 'J d k \2dk k fc 4 " ) (q- kj 



1 



7r 2 -\nq 2 /^(\nq 2 /^ + l); (11) 



R 



d 4 fc InfcVA 2 
fc 4 (<j-fc) 2 



d 4 fc 



/l d ln 2 k 2 /jl 2 
\4dk fc 4 
1 



Info - kfjpi 1 
(q - fc) 2 



(12) 



^ /- d4fc ln 2 (( g -fc) 2 / M 2 ) = /- d4jfe /I fc lnfc 2 /£ 2 ^ ln 2 ( g -fc) 2 / M 2 



k 4 (q-k) 2 J \2dk fc 4 J (q - fc) 2 

1 

^2 



- tt 2 ^ (mV/AW//* 2 + 1) +4C(3)) ; (13) 
p ^ ln(( g -fc) 2 / M 2 )lnfc 2 //2 2 _ /• , 4 , /l d > 2 fc 2 /^ ln(g - fc) 2 //z 2 

1 ' fc^^-fc) 2 - J dk [idk P J (i - fc) 2 

= tt 2 ! (V/ M 2 (l + lng 2 /A 2 + ± In 2 q 2 //) + 2f(3) ) . (14) 

Here fi is an IR renormalization parameter and Q{z) Riemann zeta-function. 

3. To calculate the two-loop master massless diagram (see Fig. 1) strictly in four 
dimensions 



Figure 1: Two- loop master diagram. 



d 4 fcd 4 Z 



J ^ ' ' k 2 (k-q) 2 (k-l) 2 l 2 (l-q) 2 ^ 
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let us apply four-dimensional IBP 



diMH w a { ^{k- q f{k-im-^ ) = °- (16) 

Formally, this is the same identity that was used within the standard IBP method 
based on dimensional regularization. Using relations 2(fe — 1)1 = k 2 — (k — I) 2 — I 2 , 
2(k -l)(l-q) = (k- q) 2 - (k - I) 2 -{I- q) 2 , one gets 

,4,,4, (k-l) 2 -(k-q) 2 _ 

d ^ l k 2 {k-q) 2 {k-l) 2 l 2 {l-qY [L7) 

Although this integral is UV and IR convergent one obtains divergences of both kind when 
considering separately two terms in the numerator. It is possible however to separate 
them by introducing (in advance) R and R operations. For example, for the i?-operation 
one uses the following identities (with j = 1, 2, . . .): 

J Ak\n J l 2 /fi 2 Ii{k,l,q) = J dfclnfcV^ln^V^i (q^- + l^-X-4\ll(k,l,q), 

(18) 

where II is a homogenous function of k, I, q of degree A. Then one obtains 

R *J l d4kdH kHk- q mi- q r =R *J S dAkdH k 2 {k-i) 2 i 2 (i- q y (19) 

where R* = RR. 

Let us now apply four-dimensional IBP with additional logarithm: 
, 4 , , 4 , d f [k- l) a \nl 2 /ft 2 



dM l dr a {k 2 (k-q) 2 (k-l) 2 l 2 (l-q) 2 ?=°- (20) 



After that one comes to 



,4, i 5 /\ 4 >((5-07m 2 ) 



J(q) = R / d 4 fc— — R / d 4 l 



k 2 {q-k) 2 J l 2 (q-l) 4 

-r R 1 ,1V,,,,, 1 * faH^r-R* I f^an- ln ^ 2 



k 2 (q-k) 2 J l 2 {q~l) 4 J J k 2 {k-l) 2 l 2 {q-iy 

R* f ^4,^4, MG-g) 2 /M 2 ) r91 s 

J J l k 2 (k-i mq -ir (21) 

For evaluation of recursively one-loop Feynman integrals involved, one uses (|6|,^0|-|L4]) 
and arrives at the well-known result |]l3|] 

J(q) = 6C(3)n 4 /q 2 . (22) 

4. As in the above example of the master two-loop diagram one can apply four- 
dimensional IBP with differential R and -/^-operations at least in every situation where 
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IBP within dimensional regularization is applied. To do this one uses almost the same 
identities as in the case of dimensional regularization. Roughly speaking, negative powers 
of e = (4 — d) /2 (where d is the space-time dimension) are replaced by some powers of 
logarithms of the loop momenta. To be able to consider separately individual terms in 
these identities one removes artificial UV and IR divergences in advance (in the case 
of dimensional regularization they were automatically regularized). To perform further 
calculations one needs a table of one-loop integrals with logarithms (in the considered 



Note that within analytic regularization one fails to use IBP for calculation of the 
master diagram, as it is possible within dimensional regularization and, strictly in four 
dimensions, in the presented technique. This can be explained by impossibility of cancel- 
lations of the type p 2 jp 2 — 1 in the integrands (because one obtains p 2 in denominators in 
a power dependent on the regularization parameter). Therefore, although the table of the 
one-loop integrals was derived with the help of analytic regularization, the logarithms 
involved do not spoil these cancellations (graphically, corresponding to contraction of 
lines to a point) which lead to recursively one-loop integrals in the right-hand side of 



This fact looks promising for construction of differential renormalization schemes 
compatible with the gauge invariance because Ward identities are finally connected with 
the possibility of the above cancellations in renormalized Feynman integrals. It turns out 
that it is possible to modify the prescriptions of differential renormalization of refs. |) 
in such a way that the gauge invariance will be preserved automatically (i.e. without 
adjusting finite counterterms) at least in the Abelian case.^j This assertion will be justified 
in a separate paper. 

Finally, let us stress that for applications in theories with chiral and super symmetries, 
the four-dimensional technique looks more preferable. 
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